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I. INTRODUCTION 



The limit of infinite spatial dimensions has turned out to be a natural starting point for 
obtaining sensible approximate |I|J^J^,|], and even essentially exact solutions of models 
of highly correlated electronic systems. In this limit the dynamics of the system become 
essentially local |^ which considerably simplifies the task of calculating quantities of interest 

mm- 

In the present paper we want to extend our previous study of the Hubbard Hamiltonian 

H = ^ij {'^la^j^a + h.C.) - flj^nia (l) 

{ij),a iff 

i 

in the limit of infinite spatial dimensions d = oo. The notation in (1) is the standard one 
and the limit d oo has to be taken such that t*^ = (i(tk)k = constant. In one of us 
demonstrated that the one-particle Green's function, or equivalently the proper one-particle 
self energy of the model (1), in this limit is obtained from the equation 

G,,iz) = I duA,{u) ^ — - = G{z) . (2) 

J Z — U) — 6 — 1-iyZj 

Here, the Green's function Q{z) is the solution of a single-impurity Anderson model with an 
effective hybridization given by [0 

A{z) = -^ + n^)~z-ii (3) 

and Ao(e) denotes the free one-particle density of states (DOS). Note that for a given site i 
equation (3) defines an effective potential due to the presence of the lattice. Equations (2) 
and (3) thus constitute the "natural" mean field theory for the Hubbard model (1) 0]. 

This mean field theory is of course independent of the lattice structure. For reasons 
of convenience, however, we shall concentrate on a simple hypercubic lattice with sites 

oo d 

and transfer along the ci-coordinate axes only, i.e. tk = — 2 ^ X! cos(m/c„). The latter 

m=l n=l 
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assumption obviously oversimplifies the situation when one wants to consider transfer beyond 
nearest neighbours, but it has the advantage that the free single-particle DOS 



when t*^ = J2 t*m^ = 1- The latter convention will set the energy scale used for the remainder 
of this paper. For nearest-neighbour transfer and q = (vr, tt, . . .) one has perfect nesting. 
However, any t2 7^ destroys this property and thus allows us to continuously bias quantities 
which depend on the perfect nesting like magnetic instabilities. 

The situation with nearest-neighbour transfer only was explored in refs. 5&7 using a 
quantum Monte Carlo (QMC) method to solve the impurity Anderson model. We could thus 
obtain essentially exact results for the model (1) and discuss magnetic and single-particle 
properties for a variety of model parameters and temperatures. The results at half filling 

= 1 can most conveniently be presented in the phase diagram in Fig. 1: For small val- 
ues of U and high temperatures one finds a paramagnetic metal with correlation enhanced 
Fermi-liquid parameters. By increasing U for a fixed temperature a crossover through a 
semi-metallic like (shaded region) into a Mott-Hubbard like phase with exponentially re- 
duced DOS at yU takes place. Note, that one never finds a true gap in the DOS for this 
"phase" . Nevertheless, transport and thermodynamic properties will essentially behave like 
an insulator. By lowering the temperature for fixed U , one encounters an antiferromagnetic 
transition which is connected with a gap in the one-particle DOS due to the cell-doubling 
associated with the antiferromagnetic state. As mentioned earlier, this phase can be shifted 
to lower temperatures or even be completely suppressed by magnetically frustrating the 
system by a finite ^2- In this case the dotted line in the antiferromagnetic region in Fig. 1 
becomes important. It visualizes the behaviour of the "MI" -crossover for a fixed U when 
the temperature is lowered and shows an interesting and unexpected reentrance behaviour. 




(4) 



acquires the simple Gaussian form 0] 




(5) 
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As we will discuss later, this is connected to a competition between the "Mott-Hubbard" 
phase and the Kondo effect also present in this model [0. 

This Kondo effect is apparent as the temperature approaches some small energy scale. 
There are two ways we may define such an energy scale. First, we can obtain an "effective 
impurity" Kondo-scale Tk from the self-consistently embedded impurity by calculating the 
screened local moment on the site i, Tx{T)ii, and then extracting Tk by fitting this result to 



the universal numerical renormalization group results of Krishnamurthy et al. [III]. We find 
that Tk defined in this way is strongly temperature dependent, increasing with decreasing T 
0. This temperature dependence of Tk can easily be understood - and is in fact expected 
- from the obvious temperature dependence of the effective medium defined by (3). It also 
explains the observed rather fast disappearance of the Abrikosov-Suhl resonance (ASR) in 
the one-particle DOS, associated with this Kondo screening, with increasing temperature 
and its total absence in the semi-metallic region of the phase diagram. Second, we may 
also identify a scale To as where the ASR in the density of states reaches half its maximum 
value. This latter energy scale appears to be more physically meaningful since it shows up 
in physical quantities like specific heat and resistivity. 

Away from half filling the Mott-Hubbard phase is immediately replaced by a Fermi 
liquid with strongly enhanced Fermi liquid parameters. More precisely, we found a narrow 
resonance at n for low temperatures which leads to the observed enhanced quasiparticle 
mass. This resonance could again be traced to a Kondo-screening of the local moments with 
a dynamically generated low temperature energy scale - Tq - connected to it. The magnetic 
transition is also found to be suppressed upon doping. Finally, for greater than 20% doping, 
correlation effects become less important and the system basically behaves like one would 
expect from standard perturbation theory. 

The remainder of this paper is split into three parts. First, we will compare different 
approximation schemes to the QMC results. The main reason is that QMC is rather time 
intensive and becomes problematic for large values of U and inverse temperature (3. Also, 
by virtue of the method, the QMC process gives all dynamical quantities as function of 



Matsubara rather than real frequency and one has to use e.g. maximum entropy methods 
to analytically continue these results to real frequencies. Although this is straightforward 
for densities of states, it proves problematic for quantities like the one-particle self energy. 
On the other hand, several physical quantities need this real frequency dependence as input. 
As we will show, a good approximation scheme for this purpose is given by a self-consistent 
perturbation theory developed for the single- impurity Anderson model (NCA) In the 

second part of the paper we use this approximation to calculate free energy, specific heat, 
resistivity and optical conductivity for the model (1) for the paramagnetic phase. Finally, a 
discussion will close the paper. 



II. COMPARISON OF DIFFERENT METHODS 

One major problem in using the QMC approach to calculate physical quantities is the 
rather large amount of computer time one has to invest to obtain results for one particular 
set of parameters. Especially for thermodynamic properties, where one has to adjust the 
chemical potential to maintain a fixed filling, it is difficult to calculate a temperature series. 
It is thus clearly desirable to have some different methods to solve the Hamiltonian (1) or, 
equivalently, the single-impurity Anderson model. 

The most straightforward idea is to use standard perturbation theory in U. This is 



known to work rather well for the symmetric single-impurity Anderson model |]12|,|13|, and 
one thus may expect it to be a reasonable approximation at half filling and for small values 
of U. Away from the symmetric point it is known that at least the lowest order does not 
reproduce the correct occupation number ||^. Nevertheless it is a simple method and it 
is surely worthwhile to outline its region of applicability. It also has the advantage that 
it automatically fulfills Fermi-liquid sumrules. A rather complete discussion up to second 
order in U has been reported by Menge and Miiller-Hartmann [0,^. Since it has been 
pointed out by Georges and Kotliar that these results are not qualitatively much different 
from the lowest second order result with Hartree self consistency |T^,[TB|, we shall use the 
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latter approach here. 

The most successful approximate methods for dealing with highly correlated electron 
systems have been developed for the single-impurity Anderson model by choosing the mixing 
term as perturbation |T^. Unfortunately, the price one has to pay for leaving the Coulomb 
interaction intact is that the standard methods of perturbation theory fail. This problem 
can be nicely circumvented for the impurity problem, leading to well defined and controlled 
approximation schemes like the so-called NCA This approximation is known to 



work well when the physics of the system is dominated by spin fluctuations ||T8|JT9(] but 



fails when charge excitations become important. In this respect it may be viewed as an 
approach complementary to standard perturbation theory. In addition, the NCA tends to 
violate Fermi-liquid properties for temperatures much lower than the smallest energy scale 
in the problem |T^. However, the NCA is nevertheless quite reliable over a large interval 



of parameters including temperature 0. Since the solution of the model (1) foTd = oo 
essentially reduces to the solution of a single-impurity Anderson model it is natural to 
adopt the NCA for this problem. 

Another natural attempt is to extend the perturbation theory with respect to mixing 
directly onto concentrated systems. In this case, however, the missing features of standard 
perturbation theory complicate the problem considerably [^] and a controlled approxima- 
tion (like in the impurity case) presently does not exist. With the use of some ad-hoc 
assumptions it is nevertheless possible to set up an approximation for this problem, too. 
These theories are originally designed for the periodic Anderson model and known in litera- 
ture as XNCA 1^ and LNCA p^j23[| . Recently, one of the authors has shown that any such 



theory for the periodic Anderson model can be readily employed for the Hubbard model (1), 
too p^j25|] . In order to obtain an idea about the quality of these approximations we include 
the LNCA in our comparison. 

The single-particle density of states for the Hubbard model (1) at half filling n = 1 for 
several values of U at an inverse temperature (3 = 7.2 is shown in Fig. 2 for the different kind 
of approaches discussed before. Let us first outline the general features of the DOS as they 
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appear from the QMC results: In all cases one finds two prominent peaks at roughly ±C//2 
which have to be identified with charge excitations on and off the local levels. In addition 
there is a pronounced resonance at /j, for small values of U due to coherent movement of 
the particles in the system. This feature is suppressed when U is increased and eventually 
a pseudo-gap opens at /i. 

In comparing the different approximations to the QMC results the first thing to note is 
that the overall agreement between QMC and NCA, apart from small differences at is 
very good. For U = 6 we did not succeed in analytically continuing the QMC results. The 
only quantity we were able to obtain here is the position of the edges of the pseudo-gap. 
These were found to be in good agreement with those predicted by the NCA. We want 
to point out that for this value of U, as generally for values U well inside the "insulator" 
phase in Fig. 1 and (3U ^ 1, the NCA does not provide stable results but tends to produce 
spurious oscillations at the gap edges. However, general structures like the width of the 
pseudo-gap are reproduced with good accuracy. Nevertheless, these instabilities prohibit a 
thorough investigation of this surely very interesting part of the phase diagram at half filling. 
We want to emphasize that this problem is not intrinsic to the NCA, but rather must be 
attributed to numerical instabilities of the computer code used to solve the NCA equations. 
The reason is that structures in the NCA equations become very sharp in this region and 
eventually cannot be resolved on a discrete energy mesh. When this occurs we approximate 
these structures as poles, which gives rise to the mentioned numerical instabilities. Note, 
that this problem does not occur outside the "insulator" phase and off half filling. A rather 
interesting point is that "poles" in the NCA begin to develop exactly when the DOS at /i 
becomes exponentially small. This empirical observation was also confirmed by QMC for 
some characteristic points in the phase diagram and eventually used to find an estimate of 
the right hand border of the crossover region in Fig. 1. 

Apparently, perturbation theory in U generally fails to reproduce even qualitatively both 
the high and the low-energy features of the DOS. Only for small U may the quasiparticle 
band at // be regarded as a good approximation to the QMC data. The LNCA on the 
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other hand looks hke a too large value of U had been used. This may be attributed to 
the approximations involved which put a strong emphasis on local correlations and are thus 
likely to overestimate residual local interactions. It also clearly overestimates the charge 
excitation bands and shows little of the finer structure near the gap edges [0, but at least 
it reproduces the general features of the DOS qualitatively correctly and accounts for the 
existence of the pseudo-gap. 

In Fig. 3 we present some typical results off half filling, namely for yU = 1 (rig ~ 0.94) and 
/i = 0.5 (ne ~ 0.8) at two different temperatures /3 = 3.6 and (3 = 14.4. The value of the 
Coulomb repulsion is U = 4. Again, QMC and NCA are in good agreement concerning the 
high and the low energy features except for Fig. 3(d), where Anca(/^) comes out much too 
large, i.e. the NCA fails to account properly for the low energy physics. This is the principal 
failure directly related to the approximations involved in the NCA [p!9[| . 

Interestingly, the LNCA gives a much weaker temperature dependence of the DOS at 
fi, indicating that this approximation underestimates the characteristic low energy scale Tq. 
This is in accordance with the observation made earlier, namely that the LNCA tends to 
overestimate the role of the local correlations. Apart from this failure the general form of 
the spectra agrees at least qualitatively with the exact result. To obtain reasonable results 
from perturbation theory, we found it necessary to fix the occupancy to the QMC value by 
adjusting the chemical potential. This given, the perturbation theory apparently becomes 
better with increasing hole concentration. It nevertheless produces features which are too 
broad and rather poor imitations of the charge excitation peaks. 

A first conclusion one may draw from these considerations is that the NCA reproduces 
most of the general features of the single particle DOS with good accuracy. However, 
Fig. 3(d) clearly shows that for some choice of parameter values the most important region 
at fi is approximated very poorly. In order to achieve a better classification of the portion 
of the parameter space where the NCA constitutes a reliable approximation to the problem 
let us substantiate the differences between QMC and NCA by looking at the quasiparticle 
weight defined by 
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where = ttT is the lowest Matsubara frequency [^]. Figure 4 displays this function for 



H = 1 and f/ = 4 as obtained from QMC (circles) and NCA (squares). Note that for these 
parameter values we expect Tq ^ 1/8 i.e. C~^(^ ^ 0) ~ 8. The agreement between 
QMC and NCA is satisfactorily, especially for temperatures T > To/2. For temperatures 
T ^ To/2 the values obtained by the NCA become too large although the order of magnitude 
is still good. Things become worse as soon as T ^ To, where the NCA produces again an 
upturn instead of a saturation. Both, the slightly too large values as well as the failure for 
T < To/5, must be attributed to the well known pathology of the NCA ||l9l. In our case, the 
most important aspect of this pathology is its tendency to give a slightly too small absolute 
value for the self energy near /i 0. While for the Anderson impurity model this behaviour 
is not important for temperatures T ^ O.lTx 0, the self-consistency process involved here 
naturally accumulates this deficiency, leading to the observed small discrepancies between 
QMC and NCA in the spectra and C^^(T). Eventually the NCA breaks down completely 
for temperatures small compared to To. However, from our comparison one can conclude 
that the results produced by the NCA are reliable for T/Tq > 1/5. 

For small U and/or far away from half filling the latter restriction makes the NCA 
obviously rather useless, because the low temperature energy scale is usually of the order 
0(1) here (see, for example. Fig. 3d with fi = 1.5 and (3 = 14.4). However, in the interesting 
region of large U and close to half filling the low temperature scale To is much smaller. In 
these cases the NCA provides a fast and consistent way to obtain information that is hard 
to access by other methods. For example, let us discuss the one-particle self energies and 
the reentrance behaviour found in the phase diagram at half filling (see Fig. 1). 

The imaginary part of the one-particle self energy is a quantity interesting in its own right 
since it provides valuable information about the low temperature behaviour of the system. 
For a normal Fermi liquid when T — 0, one expects — ImE(co' + i6) to exhibit a parabolic 
minimum at fi with a curvature and temperature dependence that is characteristic of the 



effective mass of the quasiparticles in the system. A way to obtain this latter information has 
already been discussed with the definition of the quasiparticle weight equation (8). Figure 
5 gives an impression of how the self energy behaves for some parameter values, namely at 
half filling (rig = 1) for two values of f/ = 2, 4 for a fixed temperature (3 = 7.2 in Fig. 5(a) 
and off half filling (ng ~ 0.94), for a fixed U = A and two characteristic temperatures 
/3 = 3.6, 28.8 in Fig. 5(b). While for [/ = 2 (full curve in Fig. 5a) and off half filling 
(Fig. 5b) |ImS| obviously develops a nice parabolic minimum at fi, the behaviour for f/ = 4 
at half filling (broken curve in Fig. 5a) is completely different. Here a sharp peak at n 
appears separated by a (pseudo-) gap from the continuum of particle-hole excitations. From 
general arguments [|14| it follows that — ImE(/i + i6) ~ 1/A{fi) in this case. It is clear that 



one will never obtain a Fermi liquid with this type of self energy [^]. Physically, this peak 



corresponds to an effective resonant scattering provided by the medium surrounding a given 
particle, thus localizing it by forming a bound state. It is nevertheless surprising to find such 
a structure when general phase space arguments rather suggest that particle-hole scattering 
near /i has to vanish Thus an important question is whether this structure is stable or 
may be replaced by the usual minimum for T —* 0. This leads us directly to the reentrance 
behaviour seen in Fig. 1. 

To study this interesting behaviour more closely we fix the Coulomb parameter at f/ = 3.5 
and scan the temperature from above the Ml-crossover region (T = 0.32) down to T = 0.002. 
Obviously, such a low temperature cannot be reached with QMC for this value of U. The 
results are shown in Fig. 6. One nicely sees the opening of the pseudo-gap as the temperature 
is lowered. Eventually, this pseudo-gap is destroyed by a very narrow resonance at fi which 
also signals the onset of Fermi-liquid behaviour. From the value of A{^) we extrapolate to 
a low temperature scale To ^ 1/400. Thus, the Fermi liquid that eventually emerges has 
extremely large Fermi-liquid parameters. Another question is why such an Abrikosov-Suhl 
resonance can built up from an insulator at all? Here we must keep in mind that we merely 
observe a pseudo-gap, i.e. the DOS around /i is never exactly zero and consequently will 
lead to a small but finite low temperature scale Tq. Whether the Fermi-liquid phase will win 
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depends entirely on the balance between the energy gain due to the delocalization of the 
particles in the narrow band at fi and the loss in correlation energy for the same reasons. It 
definitely seems more favourable for 3 < f/ < 4 but we cannot decide from the data available 
whether the transition line will finally intercept the abscissa at a f/ < oo or not. We must 
stress at this point that this whole scenario is valid if and only if we have sufficient magnetic 
frustration to suppress or destroy the antiferromagnetic transition appearing in the phase 
diagram. 



III. THERMODYNAMIC AND TRANSPORT PROPERTIES 

The peculiar features of the single-particle DOS and self energy discussed in the previ- 
ous section motivate a closer inspection of thermodynamic and transport properties of the 
Hubbard model (1) in the paramagnetic phase. Except for the one-dimensional model, a 
thorough study of these quantities in the thermodynamic limit was not possible yet. Pre- 
vious results from QMC simulations are usually restricted to relatively small values of 



U or comparatively high temperatures. Since they are carried out on a finite lattice with 
a discrete energy spectrum, they probably will also miss the Kondo effect if it persists in 
three dimensions. The simplifications arising in ci = oo, however, make it possible to give 
closed expressions for several quantities including the free energy, internal energy and optical 
conductivity which involve only the one-particle propagators in a simple way. We are thus 
in principle in a position to calculate these quantities exactly or, since we shall use the NCA 
to solve equations (2)-(3), obtain at least a very good approximation for them. 

Although the derivation of the expressions for those quantities is straight forward we will 
just state the final results and leave the mathematics to the appendix. To start with, the 
thermodynamic potential ^{T) is given by 

pn = Npn,^,~Y.T^T\n(^^^ . (7) 

Here, Qimp is the local free energy contribution from the effective Anderson impurity problem. 
Although, in principle, the knowledge of ^{T) provides everything one needs, it is helpful to 

11 



have an independent expression for the internal energy E{T), too. The main reason is that 
thermodynamic quantities are usually obtained by differentiating f2(T), which is a rather 
unpleasant task from a numerical point of view. In particular, the specific heat is a second 
derivative of Q, but it is a first of E{T). An expression for E{T) is given by 

EiT) = \T.j d^fi^) (^k + u;) A^^iu) + l/xiV, . (8) 

The last quantity we want to study in this paper is the conductivity. We restrict ourselves 
to the q = component, because without coupling to elastic degrees of freedom we do not 
expect the model (1) to exhibit any incommensurate charge density instability, i.e. the q = 
component will be the most important one. In this particular case, the limit d oo provides 
us with an extreme simplification, namely one can easily show (see e.g. ref. 28 and appendix) 
that the expression for aioj, T) reduces to 

a{uj) =1, j duj' j deA^{e)A{e, u')A{e, J + ^^tMl^:lM±^ . (g) 



Note that for ^ this is very similar to the result of Schweitzer and Czycholl except 
that our result is written as energy integrals and thus avoids their explicit sum on lattice 
sites that is impossible to evaluate in — > oo. 

Let us begin with a discussion of the properties of the model (1) at half filling. Due to 
numerical difficulties the NCA is currently not able to provide stable enough results in the 
interesting region just above the "MI" -crossover line. We therefore have to concentrate on 
a value just below the critical one and we found it to be a convenient choice to use f/ = 3. 
As it turns out, the behaviour found here is already close to what one may expect in the 
"insulating" region. 

Before we turn to the actual thermodynamic properties, we first want to give with Fig. 7 
an impression of the variation of the one-particle spectra with temperature. It is clear 
that the dip in the DOS at /i for higher temperatures is a poor replacement for the actual 
exponentially small DOS at larger values of U . However, together with the Abrikosov-Suhl 
resonance at low temperatures it gives a fairly good picture of the general temperature 
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dependence of the DOS even for U > Uc- Prom it we may anticipate the behaviour of the 
various thermodynamic and transport properties: Starting from high temperatures, one will 
encounter a temperature regime (e.g. t* > T ^ Tq) where the DOS mainly consists of two 
separated bands. For the entropy for example this means that it will be rather flat with a 
value reflecting the degeneracy of the states in the lower band, i.e. S ^ \n2. At the same 
time the specific heat will decrease and become very small. If there were a true gap we 
would actually expect Cy ~ exp(— /JAoap). The resistivity, on the other hand will be large 
and increases with decreasing temperature, while the optical conductivity is governed by 
the charge excitations of energy U and shows no Drude peak. 

Figure 8 displays the different thermodynamic quantities for the parameters under con- 
sideration as a function of temperature. In addition, we include for comparison some values 
of the internal energy E(T) as obtained from QMC (circles). Again, QMC and NCA are 
found to be in good agreement. It is noteworthy, that the internal energy becomes rather flat 
at T Ri 1/5. At the same time the entropy has a saddle point with a value of In 2, and 
the specific heat becomes small, as expected. The unphysical variation of S{T) and C(T) 
found in this interval must be attributed to numerical inaccuracy. In fact, by increasing the 
precision of the results we observed that e.g. the nonmonotonic variation in S{T) is reduced 
considerably while the value seems to approach 5 ~ ln2 with good accuracy. Also, the 
internal energy appears to be much more insensitive to numerical inaccuracies than the free 
energy. 

When the temperature is further lowered, we see a decrease of S{T) again, accompanied 
by a strong increase of C{T), which eventually shows a maximum. This peak in C{T) is 
a further fingerprint of the Kondo-effect in this model. Unfortunately, a further decrease 
in temperature is not possible, since the NCA-pathology becomes important for T ^ 0.01. 
From our experience this points towards a proper low temperature energy scale of about 
Tq 1/20, which is also consistent with the position of the peak in C{T) at T 7o/3. At 
present one can only infer from the knowledge of the properties of heavy fermion materials 
that C{T) ~ T/Tq for temperatures below the maximum. 
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From these results one may easily extrapolate to the behaviour of the system for U > 3. 
In fact, we mainly expect the extent of the fiat region in E{T) and S{T) to become larger, 
namely it should last roughly until a possible crossover into the Kondo regime begins. Then 
one will find a very steep decrease of both quantities again. For values of U smaller than 
f/ = 3 on the other hand the fiat region will shrink and the slope of the decrease for lower 
temperatures will become smaller. 

The resistivity for these parameter values as a function of temperature is shown in Fig. 9. 
Consistent with the DOS and the thermodynamics we first observe a semi-metallic increase at 
high temperatures which goes through a maximum and then decreases for low temperatures. 
Since we expect the system to behave like a Fermi liquid at low temperatures, the resistivity 
should follow p(T) = po + a- (T/TqY, where po = 0, Tq ~ 1/20 and a = 0(1). The calculated 
data generally follow this law when we also allow for a small intercept po ~ ~5 ■ 10^^ (see 
inset in Fig. 9). It is obvious that this negative intercept is a pure artifact and related to 
the pathology of the NCA: As discussed earlier, the approximations in the NCA tend to 
give a too small value of lm{Q^'-^^{u + iS))'^ near p. Since the self-energy is given by the 
difference between this quantity and the effective hybridization eq. (3), one will eventually 
encounter a temperature where causality is violated and the results by the NCA become 
meaningless. For the present parameter values this happens for T ^ 1/100. It is clear, that 
this breakdown will manifest itself strongest right at the minimum of ImS(ti; — i6), which 
happens to be always exactly at p at half filling due to the particle-hole symmetry. Since the 
low-temperature resistivity on the other hand is approximately just given by ImS(— i5) ||32|| , 
this violation of causality produced by the NCA leads directly to the observed unphysical 
value of Po. 

The picture for half filling is completed by the optical conductivity in Fig. 10. As 
expected, the case U = 3 already gives an idea how the "insulator" will look like: For 
high temperatures one finds a weak vestige of the Drude peak for ^ which at first is 
suppressed when the temperature is lowered. Note, however, that we always maintain a 
finite value for cr(0) consistent with the DOS in Fig. 7. At the same time the spectral weight 
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of the charge excitation peak ai u — U increases. When we further lower the temperature, 
the situation reverses. A Drude peak at a; = builds up again (see inset in Fig. 10) and 
the spectral weight at a; = t/ is decreased. In addition, a shoulder emerges at a; ~ 1. This 
feature must be ascribed to the additional excitations from the lower Hubbard band to the 
Abrikosov-Suhl resonance at Note that the spectral weight associated with this additional 
resonance is rather small. 

Let us finish the discussion of the half filled case with a look at the typical behaviour 
of the optical conductivity inside the "insulator" region in Fig. 1. Although the numerical 
problems prohibit a discussion of thermodynamic quantities here the results prove to be 
stable enough to allow the calculation of a. Figure 11 shows our results for the DOS 
(Fig. 11(a)) and cr(a;) (Fig. 11(b)) for a value of [/ = 6 and temperatures /3 — 1.44 (full 
curves) and (3 = 28.8 (broken curves). The high temperature result represents a point in the 
semi-metallic portion of the phase diagram (Fig. 1) and still shows a small peak at a; = 
due to the thermally induced states in the gap here. This feature is however completely 
lost in the "insulating phase", where only the charge excitation peak ai uj = U survives. 
We point out that the high temperature results for [/ = 3 are indeed similar to the general 
behaviour found here even though the larger DOS at of course leads to a finite value for 
cr(0) there. 

The situation off half fiUing is studied for the case U — A and Ue — 0.97. The variation of 
the DOS with temperature for these parameter values is collected in Fig. 12. Compared to 
the temperature dependence of A{0) in Fig. 7, wc observe a slower increase here, i.e. we have 
a somewhat smaller low-temperature scale Tq. From the value of A{0;T — 1/28.8) fa 2/y^ 
we extrapolate to a To ^ 1/30 for these parameter values. Although we are in principle 
able to trace the properties of the model (1) for a fixed electron density, it turns out that 
the obtainable numerical accuracy is not sufficient to get a reasonable result for the free 
energy. However, as already mentioned for the half filled case, the internal energy E(T) is 
much more well behaved and we shall concentrate on its behaviour here. The results for 
the thermodynamics are collected in Fig. 13. The features found are actually very similar 

15 



to the ones known from heavy fermion physics, as expected: We observe a maximum in 
C(T) at approximately To/3, where Tq 1/30 was read off the half-height of the DOS in 
Fig. 12. Note that in contrast to half filling we do not have a pronounced flat region in E{T) 
or S{T) here. The values for the entropy in Fig. 13 were obtained by direct integration 
S{T) — J C{T)/TdT. Taking into account the decreasing relative precision in E{T) and 
the resulting large errors in C(T) for lower temperatures we find, as expected, an entropy 
S ^\n2 associated with the low temperature peak in C{T). 

The resistivity in this case is shown in Fig. 14. As in the picture for half filling we observe 
an increase of p{T) at high temperatures which eventually goes through a maximum. For low 
temperatures we find a power law p{T) — a ■ {T/TqY with Tq f» 30 and a — 0{1) consistent 
with Fermi liquid theory. This time we do not observe any unphysical behaviour down to 
the lowest temperatures studied. This may be related to the fact, that the slight shift of the 
minimum of ImE(a; — i5) above // sufficiently reduces the influence of the NCA-pathology 
here. 

It is noteworthy that for half filling and off half filling the position of the maximum in 
p{T) does not seem to be related to Tq in a way similar to heavy-fermion systems. In fact, 
from the position of these maxima one would rather tend to rate the systems as weakly 
correlated. It is, however, important to remember that since the DOS at p is strongly 
temperature dependent, the Kondo scale itself is a function of temperature. 

Finally we present the optical conductivity for the parameter values [/ = 4 and rig — 0.97 
in Fig. 15 for some characteristic temperatures. The general structure is similar to half 
filling except that the Drude peak for a; — > continuously develops when the temperature 
is decreased. Also, the "Kondo" -shoulder found in Fig. 10 is not visible here. Again, the 
weaker temperature dependence of a{u} — > 0) points towards a smaller Tq in this case. 
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IV. SUMMARY 



In this paper we have presented a detailed study of thermodynamic and transport prop- 
erties of the Hubbard model on a infinite-dimensional hypercubic lattice. In contrast to the 
previous study we did not take into account the antiferromagnetic ordering expected 
for this model with nearest neighbour transfer, but concentrated on the paramagnetic case. 
The importance of such a study is motivated by the fact that inclusion of transfer beyond 
nearest neighbours will magnetically frustrate the system and thus depress the ordering. Be- 
sides, the effect of an antiferromagnetic transition on thermodynamic quantities like S{T) 
and C{T) is well known once their behaviour in the paramagnetic regime is known. 

The first part of this paper was devoted to a comparison between different approaches to 
the Hubbard model. As reference point we used the essentially exact QMC method discussed 
extensively in refs. 5 and 7. We found that a good description can be achieved by using 
the NCA to approximately solve the impurity Anderson model which enters the solution of 
(1) in infinite dimensions. Using the NCA, the self-consistent set of equations can be solved 
very quickly which enabled us to present a variety of quantities as functions of temperatures 
for physically meaningful parameters. One of our main results is the interesting variation of 
the entropy and specific heat in the half filled case. We believe that this behaviour can be 
viewed as generic for the strongly correlated model. This conjecture was basically confirmed 
by the qualitative similarity between our results for half filling and n = 0.97. 

Unfortunately, the NCA-approach breaks down for too low temperatures due to an in- 
trinsic violation of Fermi-liquid properties. However, our low-temperature results strongly 
suggest that a heavy-electron liquid builds up with a unique energy scale deductable from 
the variation of the various physical quantities. One thus can in principle adopt the well 
developed phenomenology for these systems to extrapolate to a consistent low-temperature 



limit for the paramagnetic phase of the infinite dimensional Hubbard model |33 



Together with our previous study, to the extent of our knowledge this represents the 
first consistent and reliable study of dynamic and thermodynamic properties of the d > 
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1 Hubbard model in the thermodynamic hmit. Although the qualitative effects of finite 
dimensional corrections are not well understood, we believe that many of the features found 
here will basically persist in at least d = 3. 
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APPENDIX A: EXPRESSIONS FOR n{T) AND E(T) 

In the following we will provide the derivation of the expressions for the free energy (9) 
and internal energy (10) for the Hubbard model in d oo. 



Let us start with the free energy. According to Baym P,|3^ one can quite generally write 
the grand canonical potential as 

/5fi(S) = $(G) - Tr(S • G) - Tr ln((G°)~^ - S) (Al) 

where shorthand matrix notations 

and 

were introduced. The functional $ is defined via the perturbation expansion of S in terms 
of G by the property 

For d —>■ oo we know that S(z) = and its perturbation expansion involves only the 

local component of G, Gii{z). This implies that $(G) = ^{{Gu}) = N^impiGu) gg] and 
finally 



(3^1{T) = iV$,^p(G',,) - iVTr(SG,,) - iVTr In(Qi) - ^Tr In 



k 



Here, fijmp is the local free energy contribution from the effective Anderson impurity problem 
for a given S [Q. 

For the derivation of E{T) we first note, that from Vt = E — TS — one obtains for 
fixed particle number E = d{(3Q) /d(3 + fiN + f3d^/df3N^. With Q = -(3'^ In Z this leads to 
the relation E{T) = (H) + fiNe. Note that this unusual form arises because our definition of 
the Hubbard Hamiltonian (1) absorbs the term —fiN. In order to calculate the expectation 
value, let us evaluate the following commutator |l35| : 



4a[Ci^, H]= tjliclha, clcj^] + cl[Ci„, c]^Ci^]) - nY ^Ifc^ 

ia i,ld,a ija 



lai ^ja\ 



i,l,a 

= Hkin -fiN + 2U = 2H- Hk,n + fiN . 
That means for (H) 

(H) = \{H,,^ + C\a[c^a, H]) - ^/xiV, 



\ Y McLcka) + ^ Yi^lalCia, H]) - ]^fiN^ 
kcr ia 

k(T ia 

'r"^ iun \ k(T ia / 



2(3 

With the equation of motion 



(A5) 



zG,,{z) = l + G,„,j {z) (A6) 



and using translational invariance, i.e. Gui^z) = j^J2Gka{z), one arrives at 

E{T) = ^YY (t^^GU^uOn) + toOnG^^iioOn) ' 1) e*""' + ^I^N, . (A7) 

^ iuJn kcr 
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Now, 

and with / duj Aiuj{ijj) — 1 

\ Y: (^u;nGk.(^u;n) - 1) e^""' = ^ E / cLu^^i^) (^^^ - l) e'^-' 

P ■ P ■ J \liOr, — LO J 

= y du;u;/(a;)Ak^(u;) 

the final equation reads 

^(2^) = ^ E / ^^/(^) (^k + ^) ^k.(^) + ^/^A^e ■ (A8) 

APPENDIX B: OPTICAL CONDUCTIVITY 

In general, the conductivity tensor is expressed via the current-current susceptibility as 

f^feiM = -^Re \^Xn,3A^ + ^<^)} (Bl) 
= ^Re{l((i,|i,))(u; + i5)} 
= Recrfe;(6i; + %S) . 

For the simple cubic lattice under consideration, the tensor is proportional to the unit tensor. 



I.e. 



^~<^) = ]^rXmim)- (B2) 

Further, the current operator for the Hubbard model (1) is given by (e = ^ = 1) 

j = E ^kJ^ka (B3) 
kcT 

where the group velocity or the particles is defined via Vk = Vtk- This leads to 
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1 1 



,,^,^._.^...„v-„,....„.-, - (B4) 

k,k ,(T ( 

This expression has the perturbation expansion shown in Fig. 16. One important imphcation 
of the hmit (i ^ oo is, that the irreducible vertex T^iUn, ioJm] i^) in Fig. 16 has to be purely 
local 1^. This means that the k-summations in the second part of Fig. 16 can be performed 
independently and thus these vertex corrections vanish, because has a different parity 
than tk [0- This means, that only the simple bubble survives and we are left with (the 
lattice constant to be taken unity) 



_ 1 1 11 



(B5) 



for nearest-neighbour transfer. 

The problem left is to evaluate the sum 

1 



Po(e) = ^EEsm^(We-tk) 



(B6) 



k I 



For this purpose, we use the method applied by Miiller-Hartmann and study the Fourier- 
transform of (B6): 

^,{s) = J po(6)e-^cie = l^^sin2(A;,)e^ 



-2ist cos k 



dk 
2^ 



sin^ ke 



-2ist cos k 



dk 
2^ 



(B7) 



d 
2 



where Ju{z) are Bessel functions. Noting that Ao{e) = J dse [Jo(2st)]'^, we find 

1 E E sm^(fc/)^(e - ^k) = ^ (po(e) + dse-^'^ Wst)]"-' M2st)^ . (B8) 



This relation is valid for a simple cubic lattice and nearest neighbour transfer in any dimen- 
sion. For the current purpose we are only interested in the limit d oo, where one can 
achieve a further substantial simplification. Noting that 2t ~ 1/ \/d ^ 1 we may approxi- 
mate J2{2st) s^t^ for large d and find for the last term in (B8) 
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J dse-''' [Jo{2st)]'^'^ J2{2st) j dss^e-''' [Jo{2st)f (B9) 

~ de^ ■ 

Since ~ 1/d this term is negligible in the limit d — > oo, i.e. po(c) = d ■ po{e)/2 and finally 

a{iu) = -^X! / deAo{e)G{e,iu}n)G{e,iu}n + ii^) = - / dedudu'Ao{e)A{e, u})A{e, ou') ^^^^ ^ ^ 

(BIO) 

Taking the real part of the analytic continuation of equation (BIO) leads to our final result 
for the optical conductivity 

a(u;) = 7r J du' J de^o(e)A(e, u;')A{e, u' + co) ^'''^"' ~ ^J'^' ^ '^^ . (Bll) 

When we collect the missing constants we find for the unit of the conductivity 

for t* leV and a — O{ao). 
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